Short-pulse photoassociation in rubidium below the Di line 
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Photoassociation of two ultracold rubidium atoms and the subsequent formation of stable 
molecules in the singlet ground and lowest triplet states is investigated theoretically. The method 
employs laser pulses inducing transitions via excited states correlated to the 55* + 5Pi/ 2 asymptote. 
Weakly bound molecules in the singlet ground or lowest triplet state can be created by a single 
pulse while the formation of more deeply bound molecules requires a two-color pump-dump sce- 
nario. More deeply bound molecules in the singlet ground or lowest triplet state can be produced 
only if efficient mechanisms for both pump and dump steps exist. While long-range l/_R 3 -potentials 
allow for efficient photoassociation, stabilization is facilitated by the resonant spin-orbit coupling 
of the Oj states. Molecules in the singlet ground state bound by a few wavenumbers can thus be 
formed. This provides a promising first step toward ground state molecules which are ultracold in 
both translational and vibrational degrees of freedom. 

PACS numbers: 33.80.Ps,32.80.Qk,34.50.Rk 
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I. INTRODUCTION 

The formation of ultracold molecules along with 
the creation of molecular Bose-Einstein condensates 
(BEC) 0, El opens the way to a new field of re- 
search Q . Methods to directly cool the translational de- 
grees of freedom of molecules have advanced significantly 
over the last few years. However, the lowest tempera- 
tures to date have been achieved by cooling atoms and 
then assembling the atoms into molecules by applying 
an external field. To this end, magnetic as well as opti- 
cal fields have been employed. The formation of weakly 
bound alkali dimer molecules with magnetic fields using 
so-called Feshbach resonances has been particularly suc- 
cessful, see e.g. H, El 0- Optical techniques, however, 
present a very general technique due to the abundence 
and diversity of optical transitions. Furthermore, laser 
fields can be manipulated offering various possibilities for 
control. In particular, molecules in their vibronic ground 
state can be created by employing a combination of light 
fields of different colors [8j- The major drawback of opti- 
cal techniques is the fairly short lifetime of the electronic 
states involved, on the order of a few ten's of nanosec- 
onds for the alkali dimer systems. The coherence of the 
control process may then be lost because of spontaneous 
emission 0, . This flaw which is particularly severe 
in view of creating a stable molecular BEC, can be over- 
come by employing laser pulses with a duration shorter 
than the excited state lifetimes. 

The use of short laser pulses offers several advantages 
in the manipulation of ultracold systems. Obviously the 
shorter timescale can overcome the losses induced by 
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spontaneous emission. In addition, a sequence of shaped 
pulses can control the photoassociation mechanism. Two 
time-delayed pulses can be employed to realize pump- 
probe experiments: The first or pump pulse creates a 
wavepacket whose dynamics is monitored by a second 
probe pulse, varying the time delay 0,0]. The shape 
of the pulse can be manipulated, for example by intro- 
ducing a controlled frequency chirp which enforces adia- 
batic following conditions 0, 0, 0, 0] . More elabo- 
rate coherent control which shapes amplitude and phase 
can also be employed. These pulses can be found us- 
ing feedback learning loops ^tJ- The underlying idea is 
to employ information on quantum interferences to con- 
struct a pulse which is able to steer the process into the 
desired direction ^^|. However, the know-how of closed- 
loop, coherent control experiments |17| cannot be trans- 
ferred to experiments on ultracold systems in a one-to- 
one fashion. This is because the dynamical timescales of 
cold systems are simply much larger than the ultrashort 
timescale of femtoseconds. The adaptation of pulse shap- 
ing techniques known from coherent control to ultracold 
systems is a subject of current investigation 0,0. In 
addition to controlling nuclear dynamics, a femtosecond 
laser can be employed for high-resolution spectroscopy in 
ultracold systems by utilizing it as a frequency comb [21j . 
These two limits have recently been combined in a spec- 
troscopic study with coherent pulse accumulation |22| . 
The appeal of such an approach lies in the fact that with 
a femtosecond frequency comb it is possible to continu- 
ously switch between the two limits, selecting the desired 
time and frequency resolution. 

Photoassociation (PA) has originally been developed 
with continuous-wave (CW) lasers, allowing for high- 
precision spectroscopy 0, 0]- PA with short laser 
pulses does not yield spectroscopic information and re- 
quires theoretical modelling. In general, PA is defined 
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as the formation of electronically excited molecules from 
two colliding atoms by interaction with laser light. In 
a subsequent stabilization step, molecules in the singlet 
ground state or in the lowest triplet state are created 
by spontaneous or stimulated emission. Photoassocia- 
tion with short laser pulses on cold systems has sparked 
interest from both experiment 0, Ha. IT^ . |20| and the- 
ory 0, 0, |25l |2|3, • Specifically, a pump-dump 
scheme has been suggested with the idea to create a 
wavepacket on the excited state and make use of its time- 
dependence for efficient stabilization 0, 0, |27j • In or- 
der to obtain an overall high efficiency, both excitation 
(photoassociation) and deexcitation (stabilization) steps 
must be efficient. PA can be understood as a vertical 
transition at the Franck-Condon distance Rc depending 
on the frequency ojp or detuning Ap. It works best at 
long range |23l |24| . Since excitation of the atomic reso- 
nance needs to be avoided, this requires a rather narrow 
bandwidth, i.e. pulses in the picosecond to nanosecond 
regime. The efficiency of the stabilization step depends 
on the topology of the potentials and the time-delay be- 
tween pump and dump pulses |2jfl . In PA experiments 
with CW lasers, the stabilization step is determined only 
by the potentials and the transition dipole moment, i.e. 
by the Franck-Condon factors of single excited state vi- 
brational levels. In a time-dependent process the excited 
state wavepacket is made up of a superposition of sev- 
eral vibrational levels and moves under the influence of 
the excited state potential. This movement can also be 
viewed as interferences in time of the superposition coef- 
ficients. Stabilization is most effective if the dump pulse 
interacts with the molecular wavepacket when it is lo- 
cated at short internuclear distances. For cesium it was 
found that up to 20% of the excited state wavepacket can 
be transferred to molecules in the lowest triplet state [27| • 

The present study employs the concepts developed in 
Refs. Ilii. 113 . PHI [lij for PA with chirped pulses and 
of Ref. |27| for the pump-dump scheme. In previous 
work on cesium, only excitation of the 0~(Ps/2) excited 
state and molecule formation in the lowest triplet state 
a 3 S+ were considered 27]. In order to analyze short 
pulse PA from the point of view of a prospective ex- 
periment, t he p resent studygoes beyond the two-state 
description [Tj, [Tj, [T]j [T^, |2]| • Formation of molecules 
in both the singlet ground state and the lowest 

triplet state a 3 S+ can then be tested. Rubidium has 
been chosen since experimental effort is under way for 
this species 0, |2(| . A comprehensive model of the in- 
teraction of two atoms with a laser field is developed 
and the role of all laser parameters is analyzed. In par- 
ticular, due to the bandwidth of the pulse, a number 
of vibrational levels of the electronically excited state 
are excited. These levels may belong to different elec- 
tronic states. For this reason all potentials which allow 
for transitions within the bandwidth of the pulse are in- 
cluded in the model. For the sake of brevity, the fo- 
cus is on the potentials correlated to the 5S + 5Pi/ 2 
atomic asymptote, i.e. to the D\ line. The corresponding 



model will be introduced in Sec. ^ Considering transi- 
tions into bS + hP\j2 potentials implies that the central 
frequency of the laser is red-detuned from the D\ line. 
The control knobs which can be varied in the experi- 
ment are the parameters of the laser pulse, i.e. its cen- 
tral frequency, spectral bandwidth, intensity and possi- 
bly frequency chirp. The variation of these parameters 
will guide the discussion of the mechanisms of excita- 
tion and deexcitation in Sees. II I II and IIVI respectively. 
The separation of pump and dump mechanisms is not 
only a convenience of the theoretical discussion. It is 
motivated by the fact that the two steps correspond to 
two different detection schemes in an experiment: The 
excitation of two colliding atoms into long-range, elec- 
tronically excited molecules would be measured as trap 
loss of the atomic cloud while the formation of molecules 
in the singlet ground and lowest triplet states would be 
monitored by R(esonantly) E(nhanced) M(ulti)P(hoton) 
I(onization) spectroscopy. Accurate information from ah 
initio calculations and spectroscopy for the molecular po- 
tentials and transition dipole moments is employed in the 
current study. Unfortunately, no such data is available as 
yet for the description of the spin-orbit coupling. There- 
fore three different model curves are employed, and the 
sensitivity of the proposed scheme with respect to the 
spin-orbit coupling is discussed in Sec. Finally, con- 
clusions are drawn in Sec. IVII In particular, the exper- 
imental feasibility of creating ultracold stable molecules 
with short pulse PA is discussed. 



II. MODEL 

Two colliding rubidium atoms in their ground state 
(5S + 5S) interacting with a laser field are considered. 
The laser excites the two atoms into an electronically 
excited state (5S + 5P) which may support long-range 
molecular bound levels. Four E + and four II (Hund's 
case (a)) states (singlet/triplet and gerade/ungerade) are 
correlated to the 5S + 5P asymptote with two of them 
repulsive and two attrative. When including the spin- 
orbit interaction (see e.g. [28|), the respective 0, 1 and 2 
(Hund's case (c)) potentials are obtained. Five of these 
are correlated to the 5S + 5Pi/ 2 asymptote, i.e. the D\ 
line. Out of these five, only the attractive potentials 
0+(Pi/ 2 ), lg(-Pi/2) an d 0g (-P1/2) which support bound 
molecular levels are considered. For the 0+ states, the 
spin-orbit coupling has resonant character [2^. l30l l3l| 
and leads to an avoided crossing between 0^(Pi/ 2 ) and 
0^(^3/2) a t short internuclear distance. The spin-orbit 
coupling therefore needs to be included explicitly: The 
Hund's case (c) potentials are obtained by diagonaliz- 
ing the Hamiltonian containing the Hund's case (a) po- 
tentials on the diagonal and the spin-orbit coupling on 
the off-diagonal |28j . The corresponding unitary trans- 
formation leads to off-diagonal elements of the kinetic 
energy which usually are neglected. For resonant cou- 
pling, however, these 'non-adiabatic couplings' have to 
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be taken into account: In the case of the two 0+(P 3 / 2 ) 
and 0^(Pi/ 2 ) curves, it couples the two vibrational series 
and leads to well-known perturbations in the spectra |29j| . 

Due to conservation of the gerade/ungerade symmetry 
in homonuclear dimer molecules, optical transitions are 
allowed to 0+ from the singlet X l H+ ground state and to 
l g and 0~ from the a 3 S+ lowest triplet state. For the in- 
teraction with the field, the dipole and rotating wave ap- 
proximations are assumed. The 0+(Pi/ 2 ), l ff (Pi/2) an d 
0~(P 1 / 2 ) states on one hand and the X 1 ^ and a 3 E+ 



states on the other hand are coupled by hyperfine inter- 
action. However, since the timescale associated with the 
hyperfine interaction is much larger than the timescales 
of femtosecond and picosecond pulses, the coupling be- 
tween the states can be neglected, and three separate 
Hamiltonians Hi, H2 , and H3 will be considered. 

The Hamiltonian Hi which describes transitions from 
the singlet ground state into the 0+ excited states is given 
by 



<T + V xl ^(R) 
Hi = I ^(R)E(t) t 




V A1 ^+(R)-fuv L 
V2W so 



\ 

V2W so 

V^CR) -Wso -huj L J 



(1) 



where T denotes the kinetic energy operator and Vi(R) 
the respective potential energy curves. The scalar prod- 
uct between the transition dipole moment and the po- 
larization vector of the field is denoted as ^ (T ( 7r )(R) for 
a (n) polarization. The laser field is characterized by 
its temporal shape and central frequency, E(t) exp(iuji,t). 
The R-dependence of the spin-orbit coupling Wso is not 
known. It is therefore approximated by its asymptotic 
value, Wso = const = 1/3APFS, which is given in terms 
of the fine structure splitting, APps = 237.5984 cm" 1 



for rubidium. In Section two different model curves 
for iyso(R) will be introduced, and the dependence of 
the results on the specific description of the spin-orbit 
coupling will be discussed. 

For the l g and 0~ states, the spin-orbit coupling has 
not got resonant character. In order to treat them on 
the same level of rigor as 0+, the spin-orbit interaction is 
nonetheless included explictly as in Eq. Q). The Hamil- 
tonians H2 and H3 then read 



and 
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so 



MR)£(t) 
f + Vs n ,(ft)- Wso-^i, 



(2) 



(3) 



Both l g and 0~ states are coupled to the lowest triplet 
state by a- as well as 7r-polarization. For simplicity the 
nuclear axis is chosen to be either parallel or perpen- 
dicular to the polarization axis of the laser field. For 
an ensemble of atoms in a magneto-optical trap (MOT) 
with no preferred axis, an average over all angles between 
the nuclear and polarization axes needs to be performed. 
This is beyond the scope of the current study. Further- 
more, rotational excitation is presently neglected and will 
be addressed in future work. 



The potentials Vi(R) have been obtained by match- 
ing the results of ab initio calculations to the long- 
range dispersion potentials (C3/R +)Cq/R +C%/R (see 
Ref. j3^| for details of the ab initio calculations). The co- 
efficients for the 55'+ 55' asymptote are found in Ref. 34], 
while the coefficients for the 5S+5P asymptote are taken 
from Ref. [3f|. The repulsive barrier of the (55* + 55*) 
potentials has been adjusted to give a triplet (singlet) 
scattering length of 99 ao (90 ao). The excited state 
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FIG. 1: (Color online) The three potentials correlated to the 
D\ line into which dipole allowed transitions may occur: Oj, 
l g , and 0g (black solid lines). Also shown are their counter- 
parts which are correlated to the D2 line and the correspond- 
ing Hund's case (a) potentials (red, blue and green dashed 
and dash-dotted lines). 
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FIG. 2: Vibrational wavefunctions of the 0^ (left) and 1 9 
(right) potentials with binding energies of 10.3 cm -1 and 
10.2 cm" 1 w.r.t. the D x line. 



The time-dependent Schrodinger equation, 

ift^|*(t)) = A(t)|*(t)>, 



(4) 



potentials in both Hund's case (a) and (c) representation 
are shown in Fig.Q] Note that at long range the 0~ (-P1/2) 
potential goes as 1/R 6 since the l/i? 3 -terms cancel each 
other. For the l s and 0~ states, the spin-orbit coupling 
influences mostly the long range part of the potentials 
while at short range the singlet (triplet) character is re- 
tained. For the potentials, however, resonant coupling 
between the -A 1 !!* and 6 3 II„ states is observed which is 
due to the (avoided) crossing at R » 9.5 ao (see inset 
in Fig. ^) . The 0+ vibrational eigenfunctions therefore 
have mixed singlet/triplet character over a large range of 
binding energies. The R-dependent transition dipole mo- 
ments are also taken from the ab initio calculations |32j | . 

The Hamiltonians, Eqs. Illl.''>ll . are represented on a grid 
employing a mapped grid method [301 136| . This allows 
to employ a fairly extended grid (i? ma x ~ 18500 a ) with 
a comparatively small number of grid points (iVgrfd = 
1027). Such a large grid is needed to faithfully represent 
the scattering continuum above the singlet ground and 
lowest triplet state potentials by box states |3 0, |53 • 
Diagonalization of the Hamiltonians, Eqs. (lll.'it . with 
E{t) set to zero gives the vibrational energy levels and 
wavefunctions. Fig. [21 shows an example of these eigen- 
functions of the 0+ and l g states. Some of the 0+ eigen- 
functions are strongly perturbed by the resonant cou- 
pling (shown on the left-hand side of Fig. see also 
Ref. I23) as compared to the regular vibrational wave- 
functions shown for the l g state. Note that the peak at 
R = 25 a.u. in the wavefunctions on the left-hand side 
of Fig. corresponds to the outer turning point of the 
Q i( p 3/z) potential. 



is solved with a Chebychev propagator. Coherent ef- 
fects resulting from laser pulses which overlap in time 
are not of interest in the current context. Therefore 
the excitation and deexcitation steps are treated sepa- 
rately. In the first step, an initial scattering state is ex- 
cited by a PA (pump) pulse. After a certain time delay 
a second (dump) pulse, suitably frequency shifted w.r.t. 
to the first one, transfers the excited state population 
back to the electronic ground state. Both laser pulses 
are assumed to have a Gaussian envelope and possibly a 
frequency-chirp, 

Ej (t) = | e ^fe4(Hf - = P)i)) , (5 ) 
V Jo 

The Gaussian standard deviation crp(u) is related to the 
full width at half maximum (FWHM) of the intensity 
profile Tptjj) of the pump (dump) pulse by Tp(rj) = 
2\/ln 2crp(c) ■ tp(D) denotes the time at which the field 
amplitude is maximum, Xp(d) ^ s the time chirp of the 
pump (dump) pulse. The stretch factor /p(d) gives 
the ratio between the pulse duration of the chirp to 
that of the corresponding transform-limited pulse, / = 
T chirp^ r ti _ j f or transform-limited pulses). Some 
of the parameters characterizing a chirped pulse are re- 
lated. While chirping stretches the pulse in time, it 
leaves its spectral bandwidth (FWHM) Aui invariant, 
Auj = 2\/ln 2r with T the Gaussian standard deviation 
of the Fourier transform of Ej(t) e~ iu> ^ t , 

E 3 {lo) = E(uj)e 3I 1 e i *i tz * iL (6) 



5 



(j = P,D) and T = 2/er tl . Time and frequency domain 
chirps are related to the spectral bandwidth via 

a j 2 4 ' v r 2 ' 

In the following, a pulse will be characterized by its cen- 
tral frequency Uj (or, respectively, detuning from the 
Di line, Aj), FWHM tj of the transform- limited pulse 
(implying its spectral bandwidth Aujj), stretch factor fj 
which gives the strength of the chirp, the sign of the 
chirp and the pulse energy. The latter is related to the 
maximum pulse amplitude Ej and to the duration <jj by 

Sj = Eoc^fKirr 2 E?j<jj 

with £o the electric constant, c the speed of light and 
r the radius of the laser beam (r = 300^m is assumed 
throughout this work). These parameters are chosen to 
correspond to pulses which can be generated from a fem- 
tosecond oscillator without amplification. This implies in 
particular pulses with FWHM of up to 10 ps and pulse 
energies of a few nano-Joule corresponding to peak in- 
tensities on the order of 100 kW/cm 2 . 

III. EXCITATION TO 5S + 5P 1/2 STATES 

The photoassociation of two free atoms to molecular 
levels of the 0+, l g , and 0~ states is first studied. The 
continuum of scattering states is represented by a finite 
number of box states 0, 0, 1^3 . These box states 
are normalized to one. To obtain the energy normal- 
ization which is usually employed for continuum states, 
they have to be weighted by the density of states of the 
box (see Ref. for details). The initial state of the 
calculations is chosen to be the s-wave scattering state 
with collision energy E co \\ equivalent to a temperature 
T = E co i\/ks = 105 /xK (with Aip, the Boltzmann con- 
stant). Such a state is typical for the conditions in a 
MOT. At these low temperatures the short-range part 
of the scattering states is independent of temperature. 
This means that the position of the nodes of the scat- 
tering wavefunctions approximately match the nodes of 
the last bound level. If the Franck-Condon radius corre- 
sponding to the central frequency of the pulse, Rc{ojl), 
is located within this region determined by the exten- 
sion of the last bound level, the excitation probability 
is roughly temperature- independent. In the following, 
7r-polarization of the laser will be assumed, and the ro- 
tational angular momentum of the molecules is taken to 
be J = 0. 

A. Excitation efficiency to the 0J, l s , and 0^ states 

For weak fields, the excitation efficiency is completely 
determined by the dipole matrix elements and by the 
detuning and spectral bandwidth of the laser. Fig. [3] 




detuning [ cm ] 



FIG. 3: (Color online) Dipole matrix elements (absolute val- 
ues squared) between an initial scattering state of two atoms 
with collision energy corresponding to T=105 /iK and excited 
state vibrational levels vs. binding energy equal to the pho- 
toassociation laser detuning. The scattering state is normal- 
ized to one (see text for further explanations). Different elec- 
tronic states are compared (the required polarization is de- 
noted in the legend). The dpms are similar for transitions to 
and l g for all detunings. For 0~ , they become comparable 
to those of the 0^ an d 1 9 states for detunings Ap > 1 cm -1 . 



shows the absolute value squared of the dipole matrix 

elements |(<? c |MR)I^t=io5uk)I 2 ( d P ms ) vs - binding 
energy E^f c or, respectively, the laser detuning, for three 
different ranges of detunings (exc = 0+,l g ,0~). Very 
close to the atomic resonance (Ap < 1cm -1 , upper 
panel), the dpms for transitions into 0~ are about two 
orders of magnitude smaller than those for transitions 
into 0+ and l g . This is due to the l/i? 6 -nature of - 
at long range as compared to 1/B? for 0+ and l g . It 
is noteworthy that this region of small detunings where 
PA is most efficient is unlikely to be accessible in an ex- 
periment with short pulses due to the bandwidth of the 
pulse. Such pulses contain spectral components corre- 
spondingto the atomic resonance, which can destroy the 
MOT |ljl|2£j. The pulse therefore has to be either far- 
detuned from the atomic line or the resonant spectral 
components have to be filtered out. This filtering can 
only be done with a finite spectral resolution which is of 
the order of 1 cm -1 . For larger detunings (Ap > 1cm -1 , 
Fig. |3 middle and lower panel) the dpms for all three 
potentials are of the same order of magnitude. The dif- 
ference between 1/i? 3 - and l/i? 6 -potentials lies only in 
the different density of molecular levels, with that of a 
l/i? 6 -potential being significantly smaller. Since for PA 
with pulses, more than one level is resonant within the 
spectral bandwidth, this density of states becomes im- 
portant for the excitation efficiency. The overall excita- 
tion efficiency for PA with short pulses will be smaller 
than one might initially expect due to the inaccessibility 
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triplet state is found for v =last-2, Pg last— 2 
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6.2 X 10 -7 rather 



TABLE I: Excited state population after the pulse (P C xc) com- 
pared to the absolute value squared of the dipole matrix el- 
ement (dpms) between the initial state and one exemplary 
vibrational level (the one which is resonant with the central 
frequency of the pulse) for transitions from the singlet ground 
(Qu) an d lowest triplet (1 3 , 0^) state. The last column lists 
the population of the last bound ground or lower triplet state 
level, respectively, after the pulse. The pulse energy is 4.2 nj 
in all cases, and 7r-polarization has been assumed. 



of the range of very small detunings. It has furthermore 
to be concluded from Fig. [3] that excitation by pulses is 
unlikely to be selective: Only very few detunings exist 
where the dpms are large for one potential and negligible 
for the other ones (cf. Ap ~ 7.5cm" 1 , Fig. [31 middle 
panel). In most cases, the bandwidth of the pulse will 
comprise resonances with significant dpms for more than 
one potential. 

Table Q] presents PA probabilities obtained from solv- 
ing the time-dependent Schrodinger equation, Eq. 
for all three Hamiltonians, Eqs. (lll.'jll . and for transform- 
limited pulses. Due to the normalization of the box state, 
the PA probability is given by the excited state popula- 
tion after the pulse, P oxc . Each of the detunings was 
chosen to be resonant with levels in both compared po- 
tentials. The corresponding dpms of the central level is 
given in the third column of Table [I] The PA probabil- 
ities for transitions to 0+ and l g at Ap = 4.1 cm -1 are 
comparable and reflect the dpms of the central level. At 
Ap = 8.6 cm , the PA probability for 0+ is an order 
of magnitude larger than for 0~, while the correspond- 
ing dpms differ only by a factor of 3.7. This reflects the 
different densities of states in a 1/R 3 and a 1/R 6 po- 
tential. The last column of Tabic [I] lists the probability 
of forming molecules in the last bound level of the sin- 
glet ground or lowest triplet state. For small detuning 
(Ap = 4.1 cm" 1 ), this probability is even higher than 
that of creating excited state molecules, a phenomenon 
which has already been observed for cesium [l|. For 
larger detuning (Ap = 8.6 cm -1 ), this probability is de- 
creased as compared to P exc but is still significant. The 
explanation of this behavior will be given below. 

The PA pulse creates a wavepacket composed of a 
superposition of several vibrational levels with time- 
dependent coefficients, on the excited electronic state. 
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FIG. 4: (Color online) The singlet ground state wavefunc- 
tion (absolute value squared) before (black solid line) and 
after the pulse (colored dashed and dotted lines). The ar- 
row indicates the position of the last node of the last bound 
ground state level. Different detunings and pulse bandwidths 
are compared: The spectral width of the pulse determines the 
size of the photoassociation window while the detuning fixes 
its position. 
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FIG. 5: (Color online) The projection of the excited state 
wavefunction after the pulse onto the Oj excited state vi- 
brational eigenfunctions (absolute values squared) vs. the 
binding energy of the vibrational levels: The excited state 
vibrational distributions reflect the detuning from the atomic 
line (peak position) and the spectral width (peak width) of the 



In addition the PA leaves a corresponding 'hole' in the 
ground state wavefunction. This phenomena is illus- 
trated in Figs. 01 and [3] for excitation from the singlet 
ground state to 0+. The ground state components 
before and after the pulse (t = tp ± 4crp) are com- 
pared in Fig. 0] for different detunings and pulse dura- 
tions corresponding to different spectral bandwidths. For 
Ap = 5.97 cm" 1 , the Franck-Condon point of the central 
frequency correponds to a maximum in amplitude of the 
ground state scattering wavefunction. It also coincides 
with the outermost maximum of the last bound level. 
The pulse cuts a hole around the position of this maxi- 
mum: Population transfer occurs within the photoasso- 
ciation window, defined by the range of distances given 
by the Franck-Condon points Rg (uj) for all frequencies 
within the bandwidth of the pulse [14( . The depth of the 
hole is determined by the intensity or pulse energy. The 
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width of the PA window is given by the spectral band- 
width of the pulse, cf. the blue dashed and red dotted 
lines in the upper panel of Fig.0](Awp = 2.94 cm -1 for 
Tp = 5 ps, Awp = 2.10 cm -1 for rp = 7 ps). Farther de- 
tuning of the central frequency from the atomic line shifts 
the PA window to shorter internuclear distances (Fig. 
middle and lower panels). At these shorter distances, 
the probability density of the ground state wavefunction 
is significantly smaller, and hence less population can be 
excited. 

Fig. El explains the relation between the PA detuning 
and the population of the last bound ground state level, 
f g ,last! for 0+. As reported in Tabled P g) i a st is larger 
than P cxc for Ap = 4.1 cm -1 , but smaller than P exc 
for Ap = 8.6 cm -1 . This is because at 4.1 cm -1 the 
center of the PA window is close to the position of the 
maximum of the ground state scattering wavefunction, as 
seen similarly in the upper panel. Since this position is 
almost identical with the position of the outermost max- 
imum of the last bound level, and since this outermost 
maximum contains about 85% of the probability density 
of this level, population transfer into this level via Rabi 
cycling is extremely efficient. At 8.6 cm -1 , the PA win- 
dow is moved toward the position of the last node of the 
last bound level (cf. Fig middle panel) , where less 
population is available to be transferred. In conclusion, 
the PA pulse is most efficient with respect to both P exc 
and P gj iast if the Franck-Condon point of its central fre- 
quency, Rc(ujl), corresponds to the position of the last 
maximum of the last bound level, and its spectral band- 
width is large enough such that the PA window comprises 
this last peak. 

Fig. displays the projections of the excited state 
wavepacket generated by the pulse onto the vibrational 
levels of 0+. The Gaussian peaks reflect the Gaussian 
shape of the pulse envelope. The peak position is deter- 
mined by the detuning of the laser, while the peak width 
mirrors the spectral bandwidth. The height of the peaks 
is entailed by the dpms which decreases for larger detun- 
ing from the atomic line (all pulses have the same pulse 
energy). Note the non-zero values of the projection close 
to the atomic line for rp = 5 ps and Ap = 5.97 cm -1 . 
Even though the spectral components close to Ap = 
are extremely small, the atomic resonance is excited due 
to its huge dpms. This emphasizes the importance of fil- 
tering out spectral components around the atomic line. 



B. Positive frequency chirp: Enhanced excitation 
probability 

The next step is to relax the condition of transform- 
limited Gaussian pulses and to introduce a frequency 
chirp where the central frequency of the pulse changes 
linearly with time. Chirped pulses have been introduced 
in the context of molecular 7r-pulses j^. The princi- 
ple underlying this concept is that if the pulse is much 
shorter than the vibrational period, the nuclear motion 
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FIG. 6: (Color online) Enhancement, i.e. ratio of excited 
state population after a positively chirped pulse to excited 
state population after the transform-limited pulse with the 
same spectral width Atop vs. stretch factor fp characterizing 
the strength of the chirp: A positive frequency chirp leads to 
more efficient population transfer for transitions to both 0„ 
and l g states. 



is not resolved by the pulse. The molecule can be treated 
as an effective two-level system, and chirping the pulse 
enforces adiabatic following conditions. PA with chirped 
pulses has been discussed within a Landau-Zener picture. 
It was found that a positive chirp maximizes population 
transfer by countering the slope of the excited state po- 
tential El [II. 

The calculations presented in Fig. show the ratio of 
the excited state population after a positively chirped 
pulse to the excited state population after the corre- 
sponding transform-limited pulse for transitions into 0^ 
and l g . The pulse has a detuning of Ap = 4.1 cm , 
a pulse energy of 4.2 nJ, FWHM of 10 ps and a cor- 
responding spectral bandwidth of Aup = 1.47 cm -1 . 
The difference between these calculations and those of 
Refs. 0, ^| lies in the description of the initial state: 
For the Gaussian wavepacket of Refs. 0, complete 
population inversion to bound excited state levels can be 
obtained, while for a scattering state as in Fig. ^popula- 
tion is transferred only from within the photoassociation 
window (see Ref. [l5[ for a detailed discussion). The 
larger enhancement for 0+ than for l g is explained by 
the different topology of the potentials leading to differ- 
ent values R c (uil)- For Ap = 4.1 cm -1 , the center of 
the PA window is found at Rc = 86.5 a.u. for 0+, but 
at Rc = 69.9 a.u. for l g , i.e. for 0„ it is close to a 
maximum, while for l g it is close to a node of the ground 
state wavefunction (cf. Fig. Therefore within the 
PA window, less population is available for transfer to l g 
than to 0+. A large amount of this little population is 
already excited by a transform-limited pulse of 4.2 nJ, 
thus chirping cannot increase it further. 



C. Negative frequency chirp: Shaping the excited 
state wavepacket 

When the chirp is negative, large frequencies precede 
small ones during the pulse. Consequently, the partial 
wavepackets at large distances are excited before the ones 
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FIG. 7: The best focussed wavefunctions for excitation to 
0Z (left) and l g (right) for transform-limited 10 ps pulses 
(middle) and negatively (top) and positively (bottom) chirped 
pulses. The A TiZ (left) and 3 I1 9 (right) components of the 
excited state wavefunctions are shown, i.e. the component of 
the state which is coupled to the singlet ground / lowest triplet 
state by the laser light: The resonant spin-orbit coupling in 
the 0j states leads to a large population at short distances 
(R ~ 20 au). tp denotes the maximum of the field amplitude, 
the stretch factors are fp =6.1 (left) and fp — 5.5 (right), 
A P = 4.1 cm -1 , £p = 4.2 nj. 



of the PA probability by a factor of three 141] ■ 

The resonant spin-orbit coupling in the case of the 0+- 
potentials ji^ causes a piling-up of population around 
R = 25 a.u. which corresponds to the outer turning 
point of the 0+(P 3 / 2 ) potential. It is in line with the 
shape of the eigenfunctions (cf. the left-hand sides of 
Figs. and In the case of the l g -potentials, the max- 
imum of population is found for R — 40 a.u. In a pump- 
dump scheme, a second pulse transfers these wavepackets 
back to the singlet ground or lowest triplet state with the 
goal of populating deeply bound levels. The positions of 
the maxima, R — 25 a.u. and R — 40 a.u., therefore 
have to be compared to the outer turning point of the 
bound ground state levels. Since levels with binding en- 
ergy larger than one wavenumber have their outer turn- 
ing point at distances shorter than R — 35 a.u., efficient 
population transfer into these more deeply bound levels 
can only be expected for 0+. Levels bound by less than 
1 cm -1 can be efficiently populated by a single PA pulse 
and therefore do not require a two-color scheme. 

Finally, the population -P g ,iast of the last bound levels 
of the singlet ground and lowest triplet states is roughly 
identical for pulses with the same frequency content, i.e. 
it is independent of a chirp. 



IV. DEEXCITATION TO X*E+ GROUND AND 
a 3 Et LOWER TRIPLET STATE MOLECULES 



at shorter distances. In an eigenstate-picture, levels with 
large vibrational period are excited before those with 
smaller vibrational period. The chirp can now be cho- 
sen to synchronize the excitation of partial wavepackets 
such that they will all arrive at exa ctly the same time at 
the inner classical turning point [Tj,[i3 ■ The value of the 
chirp can be estimated in terms of the vibrational period 
and the revival period of the level which is resonant with 
the central frequency of the pulse, i.e. in terms of the 
vibrational spectrum [T^ . 

Fig. shows how a negatively chirped pulse sup- 
presses the dispersion of the wavepacket while a posi- 
tively chirped pulse enhances the dispersion. A similar 
behavior has been reported for cesium in Refs. 0111. 
The focussing effect for a negative chirp compared to a 
transform-limited pulse is less pronounced than reported 
in Refs. 0, 0|. This is due to the larger detuning: 
Close to the dissociation limit, the vibrational periods 
increase exponentially and more correcting for dispersion 
is needed. Therefore active shaping of the excited state 
wavepacket by a negative chirp becomes less crucial with 
increasing detuning. 

A negative chirp compared to a transform-limited 
pulse can suppress the PA probability. In the range of 
detunings of a few wavenumbers, negative chirp either 
left the PA probability constant or reduced it by up to 
a factor of ten. The reduction in PA probability was ob- 
served for both 0+ and l g potentials. This is unlike the 
case of cesium where a negative chirp lead to an increase 



In order to obtain stable molecules in the singlet 
ground or lowest triplet states, a second or dump pulse 
may be applied to the excited state wavepacket. Some 
vibrational levels in the AT 1 E+ ground and a 3 S+ lower 
triplet states are already populated by the first pulse. 
However, these levels are very loosely bound. A dump 
pulse can be optimized to populate more deeply bound 
levels. The concepts developed in Ref. will be em- 
ployed, in particular the time-dependence of the dpms be- 
tween the excited state wavepacket and the bound ground 
/ lowest triplet state levels, and the idea of employing a 
narrow-bandwidth pulse, suitably detuned w.r.t. to the 
PA pulse, to achieve transfer into a single vibrational 
level. As was shown in Ref. |27|. for weak dump pulses 
the population transfer to the lowest triplet state is com- 
pletely determined by the dpms and the central frequency 
and spectral bandwidth of the pulse. The difference be- 
tween the current study and Ref. |27| is due to the stabi- 
lization mechanism allowing for efficient transfer to more 
deeply bound levels. In the study on cesium, the stabi- 
lization mechanism was most efficient for the wavepacket 
localized at the soft repulsive wall of the outer well of 
the 0g(Pz/2) state. Presently, for rubidium below the 
Di line, stabilization is afforded by vibrational levels of 
0« {Pi 1 2) which are in resonance with a level of 0+(P 3 / 2 ). 

Since coherent effects between the pump and the dump 
pulses are neglected, the excited state population after 
the PA pulse can be normalized to one. Populations cor- 
respond then directly to probabilities. Fig. [S] shows the 



FIG. 8: (Color online) The sum over all projections (absolute 
values squared) of the excited state wavepacket to bound sin- 
glet ground / lowest triplet state levels after excitation to Oj 
(black dotted lines), to l g (blue solid line, top), and to Og 
(green solid line, bottom): Deexcttation into bound levels is 
far more efficent from the 0^ state than from l g or Og . The 
transform- limited excitation pulses have rp = 10 ps (top), 
rp — 5 ps (bottom) and pulse energy of 4.2 nj. 



sum over the dpms between the excited state wavepacket 
ty e {t) and all bound singlet ground / lowest triplet state 

levels, £„»K^MR)l*e(*)>| 2 fOT ^ ^ ^ (* 
denotes the channel which is coupled by the laser field to 
the singlet ground or lowest triplet state in Eqs. (lll.'il) . i.e. 
A 1 E+ and 3 S+). The sum over projections is compared 
for 0+ and l g at a PA pulse detuning of Ap = 4.1 cm -1 , 
and for 0+ and 0~ at Ap = 8.6 cm -1 . The origin of 
time is set equal to the time when the amplitude of the 
PA pulse is maximum, tp = 0. The sum over projec- 
tions shows oscillations which reflect the excited state 
wavepacket dynamics. At later times, the oscillations are 
washed out due to wavepacket dispersion (the excitation 
pulses are transform- limited) . At both detunings, the 
projections are larger for 0+ than for l g and 0~, respec- 
tively, i.e. stabilization is expected to be most efficient 
for0+. 

The sum over the dpms does not reveal which levels 
in the singlet ground and lowest triplet states are ac- 
cessible by the dump pulse. Fig. therefore shows the 
dpms, |(^ / ,|)ti ff (R)|W e (t))| 2 ) as a function of time and 
binding energy, E v n , of the singlet ground / lowest triplet 
state levels v" (same parameters as in Fig. [HJ - For both 
l g and 0~ states, deexcitation can transfer population 
only into the last two to three levels of the lowest triplet 
state (Fig. [§]b, d). These levels are bound by less than 
0.1 cm -1 , and they are already populated efficiently by 
the PA pulse (cf. TableQJ. In contrast, for 0+ dpms with 
both the last three and more deeply bound levels are 
significant. The levels which are bound by 10 cm -1 to 
4 cm -1 have vibrational index v" = 109 to v" = 112 (the 
last bound level has v" = 120). The outer classical turn- 
ing points of the corresponding wavefunctions are located 
between R ~ 24 a.u. and R ~ 28 a.u. (as compared to 



FIG. 9: The projections (absolute values squared) of the ex- 
cited state wavepacket to all bound singlet ground / lowest 
triplet state levels vs. time and binding energy after excita- 
tion to 0+ (a,c), l g (b), and 0~ (d) (a,b: Atop = 4.1 cm -1 , 
c,d: Aujp = 8.6 cm" 1 ): Deexcitation from the l g and 0J 
states can populate only the last bound levels, while 0^ allows 
for transitions into more deeply bound ground state levels. 
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FIG. 10: The projections (absolute values squared) of the 
excited state wavepacket to all bound ground state levels vs. 
time and binding energy after excitation by negatively (a) 
and positively (b) chirped pulses to 0i (fp = 6.1, cf. the lhs 
of Fig. |7| for examples of the excited state wavepackets and 
Fig. for the projections of the corresponding transform- 
limited case): Deexcttation into more deeply bound ground 
state levels is possible after transform-limited as well as after 
chirped pump pulses. 



R ~ 170 a.u. for the last bound level), i.e. their outer- 
most maximum is located at about the same position as 
the peak in the excited state wavepackets in the left hand 
side of Fig. It can therefore be concluded that the res- 
onant spin-orbit coupling in the 0^ states provides the 
stabilization mechanism into more deeply bound ground 
state levels. It is obvious from Fig. and c that the vi- 
brational period of the excited state wavepacket as well 
as its compactness is determined by the detuning as ex- 
pected. 

The dpms of the 0+ excited state wavepacket after ex- 
citation with chirped PA pulses is shown in Fig. ^J] A 
negative frequency chirp of the excitation pulse leads to 
a compact wavepacket oscillating in the excitated state 
potential. These oscillations are reflected in the clearly 
separated peaks of the projections (Fig. HUk). After a 
positively chirped excitation pulse the wave packet is 
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spread out in coordinate space and the peaks of the pro- 
jections are smeared out (Fig. HUbb Since chirping does 
not affect the spectral bandwidth of the pulse, the ex- 
cited state wavepackets after chirped and after transform- 
limited pulses are composed of the same vibrational lev- 
els. Hence the ground state levels which are accessible 
by a second pulse are identical (cf. Fig. 1101 and Fig.[5Ji). 

The outcome of the calculations with a second or dump 
pulse are now presented. The goal is to transfer as much 
of the excited state wavepacket as possible into a single, 
more deeply bound ground state level. Calculations have 
therefore been done only for the 0+ states. The initial 
state of the example presented in the following is given 
by the excited state wavepacket after excitation by a PA 
pulse with Ap — 4.1 cm -1 and rp = 10 ps (cf. left- 
hand side of Fig. 0). The time-delay between PA and 
stabilization pulses which best achieves the goal can be 
read off the time-dependent dpms (cf. Figs. I§1 and 1 101 for 
example to — tp = 81.5 ps after the transform-limited 
PA pulse). The spectral bandwidth of the dump pulse 
is determined by the requirement of populating a single 
vibrational level, chosen to be v" — 111. Awp, needs 
to be smaller than the vibrational level spacing which is 
approximately 1.7 cm -1 for this level. This corresponds 
to dump pulses with FWHM t d > 8 ps. The detun- 
ing is chosen such that the energies of the excited state 
wavepacket and the target level are brought into reso- 
nance, Ac = —1.64 cm -1 for v" = 111, i.e. the dump 
pulse is blue-detuned w.r.t. the D\ line. The ground 
state population, P g = \ (g\^ (t))\ 2 , and the population of 
bound ground state levels, P bo und = J2 V " I (vl" l 2 > 
after a dump pulse of Tp, = 10 ps is shown in Fig.lllb as 
a function of pulse energy. Since the narrow-bandwidth 
dump pulses are resonant only with bound ground state 
levels and not with the ground state continuum, P g 
and Pbound are basically identical. The calculations of 
Fig. Illh . start from three different initial states, the ex- 
cited wavepacket after a 10 ps transform-limited PA pulse 
(dotted lines) as well as after the corresponding posi- 
tively and negatively chirped PA pulses (dashed and solid 
lines, respectively, cf. Figs. Eli and 110(1 . The highest 
ground state population is achieved after a transform- 
limited PA pulse. The positively chirped PA pulse pro- 
vides the least effective starting point for stabilization. 
This is due to the large wavepacket dispersion. Gener- 
ally, the achieved ground state population is much higher 
than expected from the dpms: Up to 50% of the excited 
state wavepacket can be transferred to the ground state. 
The vibrational distribution of the final ground state 
wavepacket, \{<j> 9 v ii\^ g {t))\ 2 , shown in Fig. ITTI: demon- 
strates that as intended this population ends up almost 
exclusively in a single vibrational level, namely v" = 111 
bound by 5.74 cm -1 . An example of a ground state 
wavepacket after the dump pulse is shown in Fig. Illb 
which confirms the pure eigenstate nature of ^> g (R;t). 
The enhanced population transfer by the stabilization 
pulse is attributed to a dynamical effect due to spin-orbit 
coupling: The excited state wavepacket has components 
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FIG. 11: (Color online) a) Ground state population P g 
(blue circles) and population of bound ground state levels 
Pbound (red triangles) after a transform-limited dump pulse 
with FWHM of 10 ps (solid lines) following excitation by a 
transform- limited (dotted lines), a negatively (solid lines) and 
a positively (dashed lines) chirped pump pulse, b) The final 
ground state wavepacket after a dump pulse with an energy 
of 8.5 nJ (corresponding to the black arrow in a), c) The pro- 
jection of the final ground state wavepacket onto vibrational 
levels (absolute values squared), i.e. the vibrational distribu- 
tion vs. binding energy after dump pulses with pulse energies 
of 8.5 nJ and 53 nJ (corresponding to the black and red ar- 
rows in a). Up to 50% of the excited state wavepacket can be 
transferred to a single, more deeply bound ground state level. 



on both A 1 E+ and 6 3 n„ states, but only the -A 1 !!* com- 
ponent is coupled to the ground state by the field. When 
the dump pulse acts, the A 1 12^ component is depleted by 
the pulse, but 'refilled' by the resonant spin-orbit cou- 
pling. Thus significantly more population is channeled 
to the ground state than expected from the dpms which 
only contain the ^4 1 S+ component. This effect is most 
pronounced when the wavepacket is focussed around the 
spin-orbit peak (cf. Fig. [7] lhs) and therefore works best 
after transform- limited or negatively chirped PA pulses. 

Transform-limited dump pulses with FWHM of 5 ps 
give similar results with slightly broader vibrational dis- 
tributions. Chirping the dump pulses decreases the pop- 
ulation transfer into more deeply bound levels: The chirp 
stretches the pulses such that the duration becomes much 
larger than the time in which the wavepacket stays at 
short distances. 



V. SENSITIVITY ON THE DESCRIPTION OF 
THE SPIN-ORBIT COUPLING 

The spin-orbit coupling having resonant character for 
the 0^ states provides a mechanism for the formation 
of ground state molecules bound by a few wavenum- 
bers. The extent to which the dump step depends on 
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the specific model of the spin-orbit coupling (SOC) is 
explored. Two model curves [H, |3!| to describe the R- 
dependence of the SOC are employed. At this time, no 
ab initio or spectroscopic data allowing for a more ac- 
curate description are available. The first model curve, 
Wi(R), has been obtained from ab initio calculations 
of the independent SOC of CS2 ^3 an d scaling with 
the ratio of the fine-structure splittings A£?fs of rubid- 
ium and cesium |39| . W\(R) displays a minimum at 
R 11 a.u., and at the crossing point of the poten- 
tials (at R ss 9.5 a.u.), the SOC is reduced to about 
65% of its asymptotic value. It was argued in Ref. [3jj 
that the reduction of the coupling at the crossing point 
is probably overestimated in Wi (R) . Therefore a second 
model curve, W2(R), was introduced showing a similar 
dependence on R as W\ (R) but with the coupling at the 
crossing point reduced to only 80% of the asymptotic 
value. The i?-dependence of both functions, W\{R) and 
W 2 {R), is displayed in Fig. 2 of Ref. [33. At long range, 
both curves become constant and equal the atomic value. 

In the following, the results of Sees. II I II and IIVI for 

Wso — const are compared to those obtained with 
Wso — Wi(R) and Wso = W 2 {R) for transitions via 
0+. The binding energies of the excited state levels are 
slightly shifted for W\{R) and W^iR) as compared to 
constant coupling. These shifts are negligible with re- 
spect to the bandwidth of the pulse. The excitation 
probability is determined by the dipole matrix elements 
between the initial scattering state and bound vibra- 
tional levels of the 0+ excited states. The absolute values 
squared of the dipole matrix elements (dpms) are shown 
in the upper panel of Fig. ^]for the three different SOC 
models. Except for single levels where the wavefunction 
is strongly modified by the coupling, the dpms are almost 
identical. This is to be expected because the overlap be- 
tween the scattering state and the excited state levels 
is biggest at large distances R where all three SOC are 
constant. For the levels where the dpms are visibly influ- 
enced, the value is reduced only by a factor > 0.5. The 
overall excitation probability should therefore not signif- 
icantly depend on the specific model of the SOC. This is 
confirmed by inspection of the excited state populations 
after the PA pulse, P CX c, listed in Table ITT1 

For deexcitation to more deeply bound molecular lev- 
els by the second pulse, it is important for the PA pulse 
to excite 0+ vibrational levels which are strongly per- 
turbed. The resonant or non-resonant character of the 
excited state wavefunctions depends rather sensitively on 
the model of SOC. Recalling Fig. [21 this is not surprising 
since the resonant character is caused by the coupling at 
shorter distances (the outer turning point of 0+(P3/2) is 
R w 20 for binding energies of a few cm -1 ) where W\(R) 
and W-2{R) differ from VFso = const. Perturbed and 
regular vibrational levels can be differentiated by their 
rotational constants B„> =< 1/(2juR) > as discussed 
in detail in Refs. [3(], I33. The rotational constants for 
the three different SOC are therefore compared in the 
lower panel of Fig. 1121 The peaks which are superim- 
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FIG. 12: (Color online) The absolute values squared of the 
dipole matrix elements | {yJT c l/x^ (R) | ¥ 3 r=i05/j.K ) | 2 (P°P) an d 

the rotational constants B v i = (l/(2/iR 2 )) of the 0j excited 
states (bottom) for three different models of the spin-orbit 
coupling. The peaks of the rotational constants correspond 
to levels which are strongly perturbed by the resonant cou- 
pling. The specific model of the spin-orbit coupling has al- 
most no influence on the dipole matrix elements. It does affect 
the position of the strongly perturbed levels which are cru- 
cial for deexcitation toward more deeply bound ground state 
molecules. Therefore the excitation probability is expected to 
be roughly independent of the description of the spin-orbit cou- 
pling. However, in order to assure an efficient dump step, the 
pump detuning needs to be adjusted such that resonant excited 
state levels are populated. 
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TABLE II: Excited state population after the PA pulse (P eX c) 
for transitions to 0j and three different models of the spin- 
orbit coupling, Wso. Also listed are the absolute values 
squared of the dipole matrix element (dpms) between the ini- 
tial state and the vibrational level resonant with the central 
frequency, and the population of the last bound ground state 
level after the pulse. The employed detunings are indicated 
by arrows in Fig. 1121 The pulse energy is 4.2 nj in all cases, 
and 7r-polarization has been assumed. 



posed on the smooth dependence of the B v i on binding 
energy (or detuning) indicate the strongly perturbed lev- 
els. For the initial detuning of 4.1 cm"" 1 (arrow labelled 
in Fig. ll2fl and a bandwidth of a picosecond pulse, per- 
turbed levels are excited only for constant SOC. In case 
of the i?-dependent coupling of W\{R) and W2(R), reg- 
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ular vibrational levels similar to those of the right-hand 
side of Fig. |21 are populated (see also Figs. 10 and 11 
of Ref. |23). In that case, it is expected that no deeply 
bound ground state molecules can be created by applying 
the second (dump) pulse. In order to investigate whether 
formation of more deeply bound ground state molecules 
is possible at all in case of /("-dependent SOC, the de- 
tuning of the pump pulse has been adjusted such that 
resonant levels are excited (arrows labelled 1 and 2 in 
Fig. El cf. also Table HTf . For larger detuning (7.6 cm -1 
and 9.2 cm -1 as compared to 4.1 cm -1 ), the vibrational 
spacing becomes larger. A smaller number of vibrational 
levels is then resonant within the bandwidth of the PA 
pulse (for Tp = 10 ps, ~ 10 levels as compared to ~ 20 
at Ap = 4.1 cm -1 ). On the other hand, the main rea- 
son for choosing a narrow-bandwidth pulse of Tp = 10 ps 
was to avoid excitation of the atomic resonance. This 
becomes less likely for larger detuning even if the band- 
width of the pulse is increased. Therefore, calculations 
for Tp = 10 ps are compared to Tp = 5 ps in Table ITT1 
Both pulses have the same pulse energy of 4.2 nJ. The 
larger number of vibrational levels resonant within the 
bandwidth of the pulse for 5 ps compared to 10 ps leads 
to a higher excitation probability. As in Sec. IIIII (cf. Ta- 
ble HJ, the excitation probabilities reflect the dpms. The 
probability to populate the last bound level of the ground 
state decreases with increasing detuning (cf. explanation 
in Sec- llTTA) . 

Fig. El shows the projection of the time-dependent 
excited state wavepacket onto the ground state level 
v" = 111, P„»=ih = |(<^„ Tn JMR)|* e (t))| 2 ,_for 
the different SOC models and different pump detunings 



bind 



= 5.74 cm -1 , cf. Fig. ||a and c). The upper 
panel confirms that more deeply bound ground state lev- 
els can only be populated if resonant levels are excited in 
0+: The maximum value of P v "=m is 2 orders of mag- 
nitude smaller for W\(R) than for Wso = const. How- 
ever, if the pump detuning is adjusted such as to excite 
resonant levels, a similar probability to populate more 
deeply bound levels is observed for W\ (R) and W% (R) as 
for Wso = const (cf. the maxima of P v "=m m ah three 
panels). Note that for larger detuning, the best overlap 
of the excited state wavepacket is obtained for v" = 110 
instead of v" = 111. The corresponding binding energy 
is i^bind 110 — 7.72 cm -1 . Therefore the resulting dump 
detuning to populate this level, Ad = Ap — E£ iT ^ 110 , 
is very small for Ap = 7.6 cm -1 and 9.2 cm -1 . In or- 
der to avoid excitation of the atomic resonance by the 
dump pulse, v" — 111 was chosen as target level instead 
of v" = 110. 

Fig. El re P or ts the results of time-dependent calcula- 
tions for the deexcitation step. The pump-dump delay 
is chosen to correspond to the first maxima in Fig. El 
{Wso = const with Tp = 10 ps, upper panel, W\{R) and 
W2{R) with Tp = 5 ps, middle and lower panel). The 
ground state population after the second (dump) pulse is 
reduced for W\(R) and W2(R) as compared to constant 
SOC. However, the amount of excited state population 
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FIG. 13: (Color online) Projection of the time-dependent ex- 
cited state wavepacket onto the ground state level v" = 111 
for three different models of the spin-orbit coupling. The 
pump pulses are the same as in Table IH1 with the detunings 
indicated in Fig. 1121 In the calculations with independent 
spin-orbit coupling, non-resonant excited state levels are pop- 
ulated in the upper panel while strongly perturbed levels are 
excited in the medium and bottom panel. Only after exci- 
tation of strongly perturbed levels, appreciable probability for 
transfer to more deeply bound ground state levels is obtained. 



W so = const, A p = 4. 1 cm 
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FIG. 14: (Color online) Ground state population P g after 
10 ps dump pulses resonant with v" = 111 for three different 
models of the spin-orbit coupling Wso (tp = 10 ps for Ap = 
4.1 cm -1 , T P = 5 ps for A P = 7.6 cm -1 and 9.2 cm -1 ). 
The specific model of spin-orbit coupling influences the exact 
amount of population which can be transferred to more deeply 
bound levels. However, the scheme is robust and efficient over 
a large range of pulse energies for all Wso ■ 



which can be transferred to v" = 111 of the ground state, 
easily reaches 20% for reasonable pulse energies for all 
three coupling models. Furthermore, a saturation of the 
transfer probability as a function of pulse energy is ob- 
served for independent SOC. This implies that the dump 
step is very robust with respect to intensity fluctuations 
of the field. 



VI. CONCLUSIONS 

Pump-dump photoassociation for rubidium below the 
Di line has been analyzed with an emphasis on exper- 
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imental feasibility. In particular, a setup such as in 
Refs. 0, |2(j was considered. All potentials into which 
population can be excited by a laser field were included 
in the model. Both potentials and transition dipole mo- 
ments were based on ab initio data and accurate long- 
range expansions. 

The first experiments on ultracold molecule formation 
with short laser pulses |19l l20j did not achieve the goal 
of creating molecules from atoms in a MOT. Therefore, 
both pump and dump steps must each be optimized in 
order to produce stable molecules in their singlet ground 
or lowest triplet state. For the PA step, an efficient ex- 
citation mechanism is provided by the long-range 1/R 3 - 
nature of several excited state potentials. However, the 
pump pulse is less efficient than one might initially expect 
from experiments with CW lasers since a rather large de- 
tuning from the atomic line is required for pulsed lasers. 
This is due to the pulse bandwidth and the constraint of 
not exciting the atomic resonance. The best compromise 
for pump pulses derived from femtosecond oscillators is 
then found for pulse durations (FWHM) of 5 ps to 10 ps 
and detunings of a few wavenumbers. 

The bandwidth Awp of the pump pulse leads to the 
concept of a photoassociation window which is comprised 
of the Franck-Condon points corresponding to all reso- 
nant frequencies contained in Auip [lj|- The excitation is 
optimal if the photoassociation window covers the range 
of internuclear distances of the last maximum of the last 
bound level of the singlet ground / lowest triplet state. 
The Franck-Condon radius corresponding to the cen- 
tral frequency (or detuning) of the pump pulse Rc(Ap) 
should then be close to the location P max of the last max- 
imum of fy = \ ast (R)- This is the case for pump pulse de- 
tunings of a few wavenumbers. If one is interested only 
in exciting atoms, pulses with positive frequency chirp 
perform best. For subsequent formation of more deeply 
bound molecules in their singlet ground or lowest triplet 
state, a pump-dump scheme is required where the opti- 
mal first pulse is transform-limited. If the detuning from 
the atomic line is small enough, the PA pulse transfers 
population also to last bound levels of the singlet ground 
or lowest triplet states, respectively. This has already 
been observed in the case of cesium [l4| . The creation of 
these extremely weakly bound molecules therefore does 
not require a two-color pump-dump scheme. 

The use of two pulses can create molecules in their 
singlet ground or lowest triplet state (bound by a few 
wavenumbers) provided that an efficient stabilization 
mechanism exists. For the Rb2 states correlated to the 
5S + bP\/2 asymptote, such a mechanism was identified 
for 0+ leading to molecules in the A^ 1 E+ ground state. 
A pump-dump scheme of photoassociation below the D\ 
line via the 0+ excited state provides then an efficient 
means to create ground state molecules bound by a few 
wavenumbers. In contrast, photoassociation via the l g 
and 0~ states will yield molecules in the lower triplet 
state which are extremely weakly bound. The excita- 
tion step from two atoms to the 0+ or l g and 0~ excited 



states is unlikely to be selective. The efficient stabiliza- 
tion mechanism for 0+ has been identified as resonant 
spin-orbit coupling |3lTl4ll ] . In a time-dependent process, 
resonant coupling leads to a dynamical enhancement of 
stabilization making the dump step even more efficient 
than deexcitation to Cs2 lower triplet state molecules 
from 0~(P 3 / 2 )- Depending on the exact description of 
the spin-orbit coupling, between 20% to 50% of the Rb2 
0+ excited state wavepacket can be transferred to a single 
vibrational ground state level as compared to 14% for Cs2 
0~(P 3 / 2 ) In order to populate a single ground state 
level, dump pulses should have a duration of 10 ps and be 
blue-detuned with respect to the atomic resonance. The 
creation of a coherent superposition of ground state lev- 
els is also possible. This requires a broader bandwidth, 
i.e. shorter duration of the dump pulses |27| . However, 
care must then be taken to avoid excitation of the atomic 
resonance. 

The 0^" excited states contain both regular and 
strongly perturbed levels. It was shown that the exci- 
tation of strongly perturbed levels by the pump pulse is 
essential for deexcitation into more deeply bound ground 
state levels. The position of the resonant levels in the 
vibrational spectrum of 0+ and hence the required pump 
detuning depends rather sensitively on the description of 
the spin-orbit coupling. Within our model, the position 
of these levels cannot be predicted accurately. However, 
the resonant levels are easily identified experimentally 
by perturbations in the level spacings or rotational con- 
stants, see e.g. Refs. [42|,|43|. Accurate spectroscopy of 
the 0+ states would improve both the potentials and the 
spin-orbit coupling and allow for obtaining a pump pulse 
from theory whose central frequency is at resonance with 
a strongly perturbed level. 

The current work has been confined to study molecule 
formation via excited states correlated to the 5S + 5Pi/ 2 
asymptote. Some conclusions can also be drawn with 
respect to the 5S + 5P 3 / 2 asymptote, i.e. the D 2 line. In 
that case, four attractive potentials into which transitions 
can be induced by the laser field, exist: 0~, l g , 0+ and l u . 
All four potentials scale as 1/P 3 at long range providing 
an efficient PA mechanism. 0~(P 3 / 2 ) is known from PA 
with a CW laser to provide an efficient stabilization 
mechanism due to the soft repulsive wall of the long- 
range well. However, a technical difficulty might prevent 
this route to be feasible: The depth of 0~(P 3 / 2 ) well is 
only about 28 cm , and pulse shaping over such a small 
frequency range is hampered by the spectral resolution 
of the pulse shaper. 

An alternative route to ground state molecule forma- 
tion might be provided by employing a femtosecond fre- 
quency comb [45J which can be operated in the picosec- 
ond to nanosecond regime. This would avoid the problem 
of exciting the atomic resonance due to large pulse band- 
width and allow for smaller detunings to be used. The 
pump step could thus be significantly enhanced, recover- 
ing the efficiency of PA with a CW laser. 

In the present study, pump-dump photoassociation has 
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been discussed in the frame of a two-atom picture, with 
zero angular momentum and no transfer of angular mo- 
mentum from the light to the molecule. One pair of 
pulses for excitation and stabilization was considered. 
The calculation of absolute molecule formation rates is 
beyond the scope of the present study. Such knowledge 
is important to estimate the laser intensity required to 
produce a detectable number of ground state molecules. 
In order to estimate absolute rates, two questions need 
to be addressed: (i) averaging over thermal (translational 
and rotational) and angular distributions of the iV-atom 
system, and (ii) the problem of the repetition rate. An- 
swering the first question requires the solution of the 
pump-dump dynamics for several initial scattering states 
(including higher collisional momentum than / = 0) and 
weighting the results with the Boltzmann factors and the 
angles between the molecular and polarization axes. The 
problem of the repetition rate stands for the fact that 
the second pair of pump and dump pulses will act on a 
different initial state than the first pair. In particular, 
the ground state population within the resonance win- 
dow might be decreased. It furthermore has to be en- 
sured that the second pair of pulses does not destroy the 
molecules created by the first pair. In order to address 
these questions, the dynamics between the pulse pairs 
has to be solved. These questions constitute the subject 
of ongoing work. 

In view of prospective applications, vibrationally ex- 
cited ultracold molecules are not adequate. The next 
goal is therefore to create ultracold molecules in their 
absolute rovibronic ground state. Molecules at temper- 
atures below 1 mK are obtained by assembling them 
from ultracold atoms due to interaction with an exter- 
nal field. The molecule formation process is coherent for 
magnetic 0, IS or optical 46] Feshbach resonances and 
for photoassociation with short laser pulses. Feshbach 
resonances and photoassociation with a single pulse lead 
to molecules in one of the last bound ground state levels, 
just below the dissociation limit. Such molecules typi- 



cally have bond lengths on the order of 100 clq and are 
bound by less than 0.1 cm -1 . At these tiny binding en- 
ergies, singlet and triplet character are mixed due to hy- 
pcrfine interaction. Two-color pump-dump photoassoci- 
ation allows for creating molecules which are more deeply 
bound. Bond lengths of 25 ao (15 ao) an d binding en- 
ergies of 5-10 cm -1 (110 cm -1 ) were found for rubidium 
(cesium p^h For subsequent Raman- type transitions to 
v — 0, molecules with bond lengths of 15-25 ao provide a 
much better starting point than the very extended Fesh- 
bach molecules: Due to better Franck-Condon overlaps, 
a smaller number of Raman steps and a reduced pulse en- 
ergy will be required • The advantage of the 0+ route 
as shown in this study for rubidium and of resonant spin- 
orbit coupling in general @ is that molecules in their 
singlet, i.e. their absolute electronic ground state are 
created. Therefore a two-color pump-dump scheme with 
picosecond pulses provides an efficient first step toward 
obtaining ultracold molecules in their absolute ground 
state. 
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